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PREFACE 


This  Memorandum,  prepared  for  the  Advanced  Research  Projects  Agency, 
is  part  of  a  study  of  those  phenomena  which  affect  the  performance  of 
optical  or  infrared  reconnaissance  and  guidance  equipment.  The  objec¬ 
tive  of  these  studies  is  to  provide  sufficient  understanding  for  the 
system  analyst  to  compute  performance  estimates  under  various  opera¬ 
tional  conditions. 

In  light  of  some  recent  experimental  meast"-ements  of  the  lateral 
phase  coherence,  the  results  of  RM-6266-ARPA  have  been  extended  to  in¬ 
clude  a  theoretical  treatment  of  the  phase  structure  function.  The 
agreement  between  theory  and  experiment  is  found  to  be  very  good  and 
supports  our  theory  of  the  mutual  coherence  function  (MCF) . 

A  quantitative  understanding  of  the  effect  of  atmospheric  turbu¬ 
lence  in  reducing  the  lateral  coherence  of  an  initially  coherent  wave- 
front  is  required  for  the  prediction  of  the  performance  of  various  de¬ 
vices  employing  lasers  for  target  acquisition  or  guidance  in  tactical 
missions.  The  MCF  is  of  fundamental  importance  in  these  applications 
because  it  determines  (a)  the  limiting  resolution  obtainable  along  an 
atmospheric  path,  (b)  the  amount  by  which  a  finite  beam  spreads,  and 
(c)  the  atmospherically  limited  signal- to-noise  ratio  using  heterodyne 
detection. 
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SUMMARY 


The  most  commonly  used  expression  for  the  wave  structure  and  mu¬ 
tual  coherence  function  (MCF)  for  an  optical  wave  propagating  in  a 
turbulent  atmosphere,  which  is  based  on  an  unphysical  extrapolation 
of  the  Kolmogorov  spectrum,  is  shown,  in  general,  to  be  incorrect. 
Along  an  atmospheric  path,  with  specified  turbulence  parameters,  the 
new  M'JF  is  shown  to  imply  greater  resolution,  less  beam  spreading,  and 
greater  heterodyne  signal-to-noise  ratios  than  indicated  by  previous 
calculations.  These  comparisons  are  made,  and  it  is  shown  that  the 
percentage  errors  in  the  previous  calculations  increase  with  decreas¬ 
ing  propagation  paths.  In  particular,  where  it  was  previously  thought 
that  the  atmosphere  limited  the  effective  coherent  detection  size  in 
heterodyne  detection  at  all  ranges,  the  present  calculation  reveals 
that  over  sufficiently  short  paths,  there  is  no  size  limit  imposed  by 
the  atmosphere. 
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I .  INTRODUCTION 


The  mutual  coherence  function  (MCF) ,  defined  as  the  cross-corre¬ 
lation  function  of  the  complex  fields  in  a  direction  transverse  to  the 
direction  of  propagation,  is  the  quantity  that  describes  the  loss  of 
coherence  of  an  initially  coherent  wave  propagating  in  a  turbulent  me¬ 
dium.  As  a  result,  the  mutual  coherence  function  is  important  for  a 
number  of  practical  applications.  It  determines  the  signal-to-noise 
ratio  of  an  optical  heterodyne  detector,  the  limiting  resolution  ob¬ 
tainable  along  an  atmospheric  path,  and  the  mean  intensity  distribution 
from  an  initially  coherent  wave  emanating  from  a  finite  aperture.  For 
a  medium  characterized  by  an  index  variation  that  is  a  gauss ian  random 
variable  with  zero  mean,  the  mutual  coherence  function  is  given  by 
exp  (-§D),  where  D  is  the  wave  structure  function.  ^ 

This  study  first  demonstrates  that  the  most  commonly  used  expres- 
(1-3) 

sion  for  the  wave  structure  function  (given  in  Eq.  (5)  on  p.  4) 

is,  in  general,  of  limited  validity.  This  expression  was  derived  on 

the  assumption  that  the  Kolmogorov  spectrum^  of  index  of  refraction 

fluctuations  can  be  extrapolated  to  arbitrarily  small  wave  numbers,  K, 

for  the  purpose  of  computing  the  structure  function.  The  apparent 

justification  for  using  the  Kolmogorov  power  law  for  K  <  L  ^  (where  L 

~  o  o 

is  the  outer  scale  of  turbulence)  is  that,  although  the  extrapolated 
spectrum  diverges,  the  integrals  necessary  to  compute  the  structure 
function  from  the  spectrum  remain  convergent.  However,  the  sensitivity 
of  the  result  to  the  divergent  spectrum  proves  to  be  considerable.  Al¬ 
though  Tatarski  (Ref.  1,  pp.  33-34)  has  cautioned  his  readers  against 
an  unlimited  extrapolation  of  the  Kolmogorov  spectrum,  Eq.  (5)  has  been 
based  on  just  such  an  extrapolation. 

In  Section  II  arguments  are  given  against  the  use  of  an  unbounded 
spectrum  (e.g.,  the  extrapolated  Kolmogorov  spectrum),  and  the  insensi¬ 
tivity  of  the  wave  structure  function  to  any  bounded  spectrum  is  demon¬ 
strated.  Based  on  this  insensitivity,  a  modified  von  Karman  spectrum 
(which  levels  off  for  K  <  Lq  )  is  used  to  examine  the  effect  of  the 
outer  scale  on  the  structure  function.  The  present  analysis  predicts 
a  transverse  (p)  dependence  of  the  plane  wave  structure  function  given 
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by  p^^[l  -  0.8(p/Lo)^^]  (where  i.Q  =  L.^/2-tt)  for  p  in  the  Kolmogorov 

inertial  subrange,  l  «  p  «  L  .  Varying  the  form  of  the  spectrum  for 

-1  °  °  1/3 

K  <  h  produces  small  variations  in  the  coefficient  of  the  (p/L  ) 
o  o 

term.  In  addition,  for  p  ~  0.1  h  ,  our  structure  function  has  a  bilog- 

o 

arithmic  slope  closer  to  3/2  than  5/3  and  levels  off  for  p  >'h  . 

(4)  ° 

Bouricius  and  Clifford  recently  measured  the  phase  structure 

function  at  0.6328  p  over  a  50-m  path  1.6  m  above  the  ground  for  trans- 

* 

verse  separations  varying  from  1  cm  to  2  m.  From  temperature  measure¬ 
ments,  they  also  deduced  the  index  structure  constant  and  the  outer 
scale  (the  latter  agreeing  reasonably  well  with  the  height  above  the 
ground) .  Their  results  are  in  agreement  with  the  3/2  dependence  dis¬ 
cussed  above  and  display  a  leveling  off  for  separations  of  the  order 
p  ~  Lq.  Buser^^  has  also  made  measurements  of  the  phase  structure 
function  at  0.6328  p  over  a  50-m  path  1.75  m  above  the  ground  for  p  * 

5/3 

22  cm;  he  obtains  a  slope  always  less  than  the  previously  predicted  p 
In  Section  III  it  is  shown  that  there  exist  three  distinct  propa¬ 
gation  distance  regimes  for  which  approximate  expressions  for  the  mu¬ 
tual  coherence  function  can  be  found.  These  formulae  with  the  respec¬ 
tive  ranges  of  validity  are  presented  in  the  Table  for  the  modified  von 
Karman  spectrum  of  Section  II. 

In  Section  IV  we  calculate  the  implications  of  our  expression  for 
the  MCF  with  regard  to  resolution  and  beam  spreading.  In  particular, 
we  show  the  coherence  length  (defined  as  the  transverse  separation  at 
which  the  MCF  is  equal  to  e  can  be  considerably  greater  than  the 
previously  accepted  value.  Comparison  with  this  value  predicts  an  in¬ 
crease  in  the  implied  resolution  and  a  decrease  in  the  implied  beam 
spreading. 

Finally,  in  Section  V  we  examine  the  implications  with  regard  to 

coherent  optical  detection.  Comparing  our  results  for  heterodyne  de- 

(3) 

tection  with  those  of  Fried,  we  predict  a  greater  long-term  average 
signal-to-noise  ratio.  In  particular,  where  Fried  finds  a  maximum  use¬ 
ful  receiver  diameter  for  all  ranges,  we  find  that,  for  distances  small 
compared  with  the  mean  field  decay  length  z£  (defined  in  Section  II), 
the  signal-to-noise  ratio  increases  indefinitely  with  receiver  size. 

*For  short  propagation  paths,  the  wave  and  phase  structure  function 
are  identical  (Ref.  1) . 
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II.  THE  SPECTRUM  AND  WAVE  STRUCTURE  FUNCTION 


The  analysis  is  based  on  the  expression  for  the  wave  structure 
function  for  the  case  of  a  plane  wave  incident  upon  a  homogeneous,  iso¬ 
tropic  turbulent  medium: ^ 


D(p,z) 


1  0 

8tt  “k'-  z 


[1  -  J  (kp)H  (K)  KdK 
o  n 


(1) 


where  k  is  the  optical  wave  number,  p  is  the  transverse  separation  at 

propagation  distance  z,  and  $n(K)  is  the  three-dimensional  spectral 

density  of  the  index  of  refraction  fluctuation.  Equation  (1)  may  be 

* 

written  in  the  form 


D(p,z) 


jo  J0(Kp)$n(K)  KdK 

io  VK>  KdK 


(2) 


where 


z 

c 


(3) 


can  be  shown^  to  be  the  propagation  distance  in  which  the  mean  field 
of  a  plane  (or  spherical)  wave,  (U),  decays  to  e  ^  from  its  value  at 
the  source  (the  angular  brackets  denote  an  ensemble  average). 

The  spectral  density  most  commonly  used  to  represent  atmospheric 
index  of  refraction  fluctuations  is  that  due  to  Kolmogorov^ 


VK> 


0.033C2K'U^3 

n 


w 


In  deriving  Eq.  (1)  we  have  assumed  that  $n  is  not  &  iunction  of 
propagation  distance.  The  modification  of  Eq.  (1)  to  include  an  explicit 
dependence  on  range  is  that  z$n(K)  +  $n(K,z/)  dz'  .  For  spherical  wave 

propagation  let  z(l  -  JQ(Kp)]  -*■  dz '  (1  -  T  fv^~t  /-'‘l 


J0(KpzVz)]. 
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valld  within  the  Inertial  subrange  ^  «  K  <<  *  \  where  A  ■  2irt 

o  o  o  o 

and  L  ■  2itL  are  the  inner  and  outer  scales  of  turbulence,  respec- 
o  o 

tively,  and  Cn  is  the  index  structure  constant. 

In  order  to  compute  D  from  Eq.  (1),  it  is  necessary  to  make  cer¬ 
tain  reasonable  assumptions  regarding  the  spectrum  outside  of  the  iner- 

(1)  (2  3) 

tial  subrange.  Tatarski'  7  and  others,  *  using  the  spectrum  of  Eq. 
(4)  for  all  K,  have  computed  the  wave  structure  function  given  by 

D  (p,z)  -  2.91k2C2zp5/3,  l  «  p  «  L  (5) 

o  *  n  o  o 


wl  lch  is  used  for  all  z,  the  apparent  justification  being  that  the  in¬ 
tegral  in  Eq.  (1)  converges.  However,  the  sensitivity  of  D  to  the  form 

of  f:he  spectrum  for  K  <  i  is  considerable.  This  can  be  seen  by  car- 

0  -1 

rying  out  the  integration  in  Eq.  (1)  for  <  K  <  •»,  which  gives 


D1(p,z)  -  2.91kVp5/3 


£  -  0.(>7(^)  +0(#)]’ 


for  l  «  p  «  L 
o  o 


(6) 


The  percentage  difference  is  [(D.  -  D  )/D.]  x  100  ■  67(p/£  )  [1  - 

1/3-1  ioi  o 

O.S7(p/i.  )  ]  and  equals  45  percent  for  p/L  "0.1.  This  example 

O  0 

reveals  the  sensitivity  of  the  structure  function,  in  the  Inertial  sub¬ 
range,  to  the  physically  unreasonable  extrapolation  of  Eq.  (4)  for  K  < 

Lq  .  Although  the  Integral  in  Eq.  (1)  converges,  an  extrapolation  would 
lead  to  divergent  integrals  for  both  the  energy  per  unit  volume  of  the 
fluctuations  and  the  distance  over  which  the  mean  field  decays  (i.e., 
z£) .  On  physical  grounds,  the  spectrum  must  begin  leveling  off  for  K 
corresponding  to  scales  large  compared  with  the  separations  over  which 
the  temperature  fluctuations  exhibit  appreciable  correlation,  with  a 
finite  upper  bound  as  K  0.  Further,  careful  examination  of  Eqs.  (1) 
or  (2)  reveals  the  insensitivity  of  the  integrals  to  the  spectral  den¬ 
sity  for  any  spectrun  which  remains  bounded  for  K  <  and  hence,  for 

the  purpose  of  computing  the  wave  structue  function,  we  suggest  the  use 
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of  any  bounded  spectrum  in  this  range.  The  spectrum  falls  off  verv 
rapidly  for  K  >  due  to  viscous  damping,  and  it  is  customary^  to 
use  a  gaussian  decay  in  this  region. 

An  example  of  a  spectrum  that  is  convenient  for  computational  pur¬ 
poses  is  the  modified  von  Karman  spectrum 


VK> 


,  -<k*  r 

0.033C  e 

(K2  +i-2,11/6 
o 


(7) 


-1  * 

which  implies  a  flat  spectrum  for  K  <  Lq  .  For  example,  substituting 
Eq.  (6)  into  Eq.  (3)  yields,  for  «  Lq, 


z 

c 


9?  (o.39k2cV/3)  1 


(8a) 


and  a  phase  structure  given  by 


D(p,z) 


2.91k2C2p5/3 

n 


l  «  p  «  L  (8b) 

o  o 


which  may  be  compared  with  Eq.  (6).  Figure  1  is  a  gra^'  -  ’  the  mean 

field  decay  length  z  plotted  versus  wavelength  X,  for  •'  "  100  cm  and 

2C 

typical  values  of  C  . 


The  quantity  27t/L0,  rather  than  1/Lq,  is  introduced  in  Eq.  (7) 
because  we  are  comparing  a  wave  number  with  a  length.  In  any  case, 
this  parameter  is  to  be  regarded  as  the  reciprocal  of  the  wave  number 
where  the  spectrum,  in  the  low  frequency  regime,  begins  to  deviate 
from  a  K“**'3  dependence.  As  has  been  noted,  the  resulting  functional 
dependence  on  p  is  insensitive  to  the  details  of  the  low  frequency  be¬ 
havior  of  the  spectrum,  but  depends  only  on  the  value  of  K  where  the 
spectrum  begins  to  level  off. 
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Substituting  Eq.  (7)  into  Eq.  (2)  and  assuming  p  »  4  yields 


4? 

r  (P.z)  -f 

r  c 


•  5  (Ji\5n  ("  Vu)u  du 

1  '  3  V  Jo  [u2  +  (»/io)2]11/6 


(p  >>  tQ)  (9a) 


The  modification  of  Eq.  (9a)  for  spherical  waves  is  given  by 

5/3 

\  du  uJ  (u)  \  f-z - 5 - 

JO  J 0  I  u  +  ( so />, 

°  J 


Da(p.z)  ■  ~  < 

S  ZC' 


,  /  \5/3  fl  5/3  .  , 

-  T  (j^-)  1  du  uJ  (u)  1  y~2 - - - S2-i11  /e., 

3  V V  JO  °  Jo  fu2  +  (sp/i_)2lll/6< 


(p  »  lQ)  (9b) 


In  Fig.  2,  the  quantity  (zc/4z)D,  which  is  independent  of  X,  z,  and  (T 
(and  4q  for  p  »  f.Q)  is  plotted  as  a  function  of  p/fcQ  for  both  plane 
and  spherical  waves.  Note  that 

a.  The  theory  predicts  a  power  law  dependence  closer  to  3/2  than 
5/3  for  separations  of  ~  0.1  and 

b.  The  structure  function  saturates  to  a  value  of  4z/z  for 

c 

p  »  &  . 
o 

The  power  law  dependence  described  in  (a)  has  recently  been  observed 

(4) 

by  Bourlciu8  and  Clifford. 
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The  normalized 
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I II .  APPROXIMATE  FORMULAE  FOR  THE  MCF 


The  mutual  coherence  function  Is  defined  as  the  cross-correlation 
of  the  complex  fields  In  a  direction  transverse  to  the  direction  of 
propagation:  M(p,z)  s  (tKr^z.U  (r2,z)>,  where  p  -  | -  r2|  .  To 
second  order  In  the  refractive  Index  fluctuations  in  general,  and  to 
all  orders  when  the  Index  fluctuations  are  a  gausslan  process,  M(p,z)  - 
exp  {-  ^l)(p,z)},  where  I)  is  given  by  Eqs.  (1)  or  (2). 

While  4t  Is  a  simple  numerical  calculation  to  compute  the  MCF  di¬ 
rectly  using  the  spectrum  of  Eq.  (7)  In  Eq.  (1),  it  Is  useful  to  have 

approximate  formulae  for  the  MCF  for  estimating  the  coherence  at  vari¬ 

ous  ranges. 

From  Eq.  (2)  It  follows  that  the  MCF  decreases  monotonlcally  from 

M(0,z)  -  1  to  M(«,z)  ■  c  c.  This  behavior  is  in  accord  with  the 

physical  picture  of  the  light  arriving  at  the  two  points  r^,  r2  being 
scattered  thi  ,  \gh  statistically  Independent  media  when  p  -  | r^  -  r, |  Is 
sufficiently  large  that  M(p,z)  -  <U<r1>U*(r2) >  -*■  <U(rj)>  <U*(r2>>  Hence 
for  z  «  z£  the  MCF  is  essentially  unity  for  all  values  of  p,  i.e., 

M(p,z)  -  1,  z  «  z  (1J) 

c 


Considering  the  Inertial  subrange,  substituting 
(1),  assuming  p  »  Lq,  and  expanding  to  lowest  order 


Eq.  (7)  into  Eq. 
in  (p/-fco>,  we  obtain 


Mj^p.z)  -  exp 


2^1  k2C2zo5/3 
—  k  Cnzp 


for  l  «  p  «  L 
o  o 

(ID 


In  order  that,  at  a  given  range,  all  of  the  transverse  separations  of 
Interest  lie  In  the  Inertial  subrange,  it  is  necessary  that  M(Iq,z)  s* 
1,  and  M(l<o,z)  «  1,  which  Is  essentially  the  condition 


z  «  z  «  z.  (12) 

cl 

where  z  Is  given  by  Eq.  (8),  and  z . ,  defined  by  replacing  i  by  t  In 
c  1  o  o 
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the  formula  for  z c,  is  a  distance  at  which  the  coherence  length  of  the 
field  is  of  the  order  of  the  inner  scale. 

For  ranges  greater  than  z^,  all  of  the  p's  of  Interest  are  small 
compared  with  the  inner  scale.  The  Bessel  function  in  Eq.  (1)  can  then 
be  expanded  in  powers  of  p/-fcQ  to  yield 

M(p,z)  a*  exp  1.72k^cjjio^3zp^,  z  »  z^  (13) 

It  follows  from  the  above  discussion  that,  for  z  «  z^,  the  plane 
wave  MCF  does  not  depend  on  the  inner  scale  and  can  be  written  as 

M  (p,z)  -  exp  [-  (p,z)l  e 
p  p  \o  c 

where  Dp  is  given  by  Eq.  (9a).  In  Fig.  3  we  compare  the  approximate 
expressions  of  Eqs.  (5)  and  (10)  with  Eq.  (14),  taking  z/zc  -  10  and 
plotting  the  results  as  a  function  of  p /tQ.  For  this  comparison,  the 
MCF  corresponding  to  Eq.  (5)  has  been  written  as 

Mo(p,z)  -  exp  [-  3.72(z/zc)(p/LQ)5/3]  (15) 

and  Eq.  (11)  as 


) 


(14) 


The  comparison  reveals  that  for  z  <<  z^,  the  field  from  a  plane 
or  spherical  wave  retains  its  transverse  coherence  to  a  greater  degree 
than  previously  predicted.  The  approximate  expressions  for  the  MCF 
with  respective  ranges  of  validity  are  summarized  in  the  Table  on  p.  12. 
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APPROXIMATE  EXPRESSIONS  FOR  THE  PLANE  WAVE 
MODULATION  TRANSFER  FUNCTION 


aThe  quantities  z  (the  distance  where  the  average  field  is 
down  by  e"“)  and  z ^  (the  distance  where  the  coherence  length  of 
the  field  is  of  the  order  of  the  inner  scale  of  turbulence)  are 
given  by  (0. 39k2C2i,5/  3)“1  and  (0 . 39k2c£«jj'  3)_1 ,  respectively. 
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IV.  RESOLUTION  AND  BEAM  SPREADING 


Two  Important  consequences  of  our  more  coherent  MCF  are  that,  for 

given  turbulence  parameters  along  a  uniform  atmospheric  path,  the  mean 

visibility  Is  better,  and  the  average  amount  by  which  a  finite  beam 

spreads  Is  smaller,  than  predicted  by  Fried. ^  Defining  p  as  the 

° 

transverse  separation  at  which  the  atmospheric  MCF  Is  reduced  by  e  , 

then  the  minimum  resolvable  length  at  a  distance  z  from  an  observer  Is 

well  known  to  be  ~z/kp  .  It  can  be  shown  thet  when  p  Is  small  com- 

o  o 

pared  to  the  size  of  the  transmitting  aperture,  the  angular  spread  of 

a  finite  beam  due  to  the  atmosphere  Is  ~  1/kp  . 

o 

As  shown  In  the  previous  section,  for  distances  small  compared 
5  3 

with  z.  (■  10  z  for  L  /l  *  10  ),  the  plane  wave  MCF  does  not  depend 
1  c  o  o 

on  the  Inner  scale,  and  can  be  written  as 


Mp(p,z) 


exp  [-  |Dp(p,z)]  e  F 


z  «  z^ 


(17) 


where  D  Is  given  by  Eq.  (9a).  The  modification  of  Eq.  (14)  for  spher- 
P 

leal  waves  Is  given  by 


Ma(p,z)  -  exp  [-  £Ds(p,z)1 


i^r)  • 

'  O  C' 


z  «  z. 


(18) 


where  D  Is  given  by  Eq.  (9b) .  In  Figs.  4a  and  4b  we  plot  p  /i.  versus 

z/z  for  plane  and  spherical  waves,  obtained  by  Inverting  F(p  /-fc  ,  z/z  ) 
c  ^  o  o  c 

e-1  and  F  (p _/■*>»  z/z  )  ■  e  ,  respectively.  On  the  same  graph  In  each 

8  0  0  C 

figure,  we  compare  our  results  with  the  MCF  obtained  from  the  structure 
function  of  Eq.  (5)  and  the  corresponding  expression  for  spherical  waves 
(obtained  by  multiplying  the  exponent  In  Eq.  (5)  by  3/8).  For  the  pur¬ 
poses  of  this  computation,  we  have  written  Eq.  (5)  In  the  form 

yp,.>  -  7.*<f-Xr)5/3 

c  o 
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The  percentage  error  in  each  case,  defined  as  [(pq  -  Pq)/pq1  x  100, 

where  p'  satisfies  M  (p',z)  ■  e  \  is  indicated  on  the  scale  at  the 
o  o  o 

right . 

As  the  plane  or  spherical  wave  source  is  approached,  the  present 

analysis  indicates  that  the  transverse  coherence  of  the  field  increases 

at  a  greater  rate  than  previously  predicted.  The  errors  increase  from 
4 

~  10  percent  at  10  z  to  100  percent  at  0.5z  ,  whereas  for  distances 

C  -1  C  5 

<  0.5z  ,  the  new  MCF  is  never  down  to  e  .  For  z  >  z.  (■  10  z  for 
c  3  i  c 

LQ/io  ■  10  ),  the  use  of  Mq  again  gives  a  poor  approximation  to  the  new 

MCF,  which  in  this  region  is  given  by  Eq.  (12). 
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V.  SIGNAL-TO -NOISE  RATIO  IN  HETERODYNE  DETECTION 


As  an  additional  Illustration  of  the  Implications  of  our  more  co¬ 
herent  MCF,  we  use  It  to  compute  the  slgnal-to -noise  ratio  In  a  coher¬ 
ent  detection  system  and  compare  the  results  with  those  computed  by 

(3)  (3) 

Fried.  The  slgnal-to-nolse  ratio  is  given  by 

1 

dx  x  K  (x)M(Dx,z)  (19) 

a  o 


where  Ag  Is  the  signal  amplitude,  n/e  Is  the  quantum  efficiency  measured 
In  electrons  per  unit  energy,  D  Is  the  diameter  of  the  collecting  aper¬ 
ture,  and 


Ko(x) 


/  \  2.1/2 
(x)  -  x(l  -  x  ) 


Then  Eq.  (19)  can  be  written 


for  0  <.z  «  z^,  where  £  »  D /lQ,  and 


(20) 


(21) 


(22) 


where  F  Is  given  by  Eq.  (14) .  The  function  ip  contains  the  dependence 
of  slgnal-to-nolse  ratio  on  the  collector  diameter  in  the  presence  of 
atmospheric  distortion,  and  Is  plotted  as  a  function  of  £  in  Figs.  5a 
and  5b  for  various  values  of  z/zc>  The  reduced  slgnal-to-nolse  ratio 
ip,  as  computed  from  Eq.  (22)  (solid  curve),  is  compared  with  the  results 
of  Fried, who  used  the  MCF  of  Eq.  (15)  (dashed  curve).  In  general, 
the  present  analysis  predicts  a  somewhat  larger  slgnal-to-nolse  ratio. 


Reduced  signal -to-noise  ratio 
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The  difference  la  more  pronounced  for  distances  small  compared  with  z  , 
where  Cs)/N  as  derived  here  increases  indefinitely  with  aperture  size. 
This  result  is  in  contrast  with  that  of  Fried,  who  derived  an  effective 
limiting  diameter  beyond  which  increasing  the  diameter  results  in  very 
little  improvement  in  (S)/N  for  all  ranges.  For  ranges  very  large  com¬ 
pared  with  z£  (but  still  small  compared  with  z^) ,  the  MCF  obtained  from 
Eq.  (1)  approaches  that  obtained  from  Eq.  (5),  and  the  difference  be¬ 
tween  the  two  signal-to-noise  ratios  tends  to  zero.  This  trend  is  il¬ 
lustrated  in  Fig.  6,  in  which  the  reduced  signal-to-nolse  ratio  is  plot¬ 
ted  as  a  function  of  z/z  for  £  ■  1. 

c 


Fig.  6 — Reduced  signal-to-noise  ratio  at  10.6  /*  in  a  heterodyne  detection 

system  a«  a  function  of  z/z  for  the  diameter  of  the  collector  equal  to  t 

c  ® 


Percentage  error 
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